Transport, atom blockade and output coupling in a Tonks-Girardeau gas 
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Recent landmark experiments have demonstrated how quantum mechanical impurities can be cre- 
ated within strongly correlated quantum gases and used to probe the coherence properties of these 
systems [JJ. Here we present a phenomenological model to simulate such an output coupler for a 
Tonks-Girardeau gas that shows qualitative agreement with the experimental results for atom trans- 
port and output coupling. Our model allows us to explore non-equilibrium transport phenomena 
in ultra-cold quantum gases and leads us to predict a regime of atom blockade, where the impurity 
component becomes localised in the parent cloud despite the presence of gravity. We show that this 
provides a stable mixed-species quantum gas in the strongly correlated limit. 



Ultra-cold gases provide an extremely versatile re- 
source of quantum matter. The possibility of applying 
external electromagnetic potentials means that one can 
create macroscopic harmonic traps and microscopic peri- 
odic structures that are sensitive to the quantum states 
of the atomic species [2J. One recent breakthrough has 
been the realisation that quasi one-dimensional quantum 
degenerate gases can be created by using strong trans- 
verse trapping potentials [3]. Such settings allow the ra- 
dial (transverse) degrees of freedom to be frozen out and 
in addition to tunc the effective one-dimensional interac- 
tion strength [3] . The resulting quantum many-body sys- 
tem of interacting bosons can be described by the Lieb- 
Liniger model, which possesses exact solutions [5j. In the 
limit of strong repulsive interactions these can be simpli- 
fied due to the existence of a mapping theorem to an 
equivalent gas of non- interacting fermions 6 . This so- 
called Tonks-Girardeau (TG) limit was recently experi- 
mentally reached independently in two separate labora- 
tories [8] and the initial experiments were followed by 
several spectacular experimental studies. These included 
the exploration of the relationship between integrability 
and thermalization [9] and also the creation and detec- 
tion of meta-stable excited states [TU] . 

Very recently Palzer et al. PQ demonstrated an out- 
put coupler for an optically-trapped quantum gas in the 
TG regime. An array of one-dimensional clouds was 
created in a two-dimensional optical lattice and then 
probed using a radio-frequency pulse to locally popu- 
late an untrapped hyperfine level. The transport prop- 
erties of this untrapped 'impurity' as it fell under grav- 
ity, passing through the parent cloud, were subsequently 
observed. This experiment constitutes a genuine open 
quantum system and also one of the first experiments 
that explores the important topic of quantum transport 
in a clean ultra-cold environment (see also [TT] ) . Here we 
propose a theoretical framework that can describe the ef- 
fects experimentally observed and predict the existence 
of a regime of atom blockade, in which the interaction 



between two components is strong enough to trap and 
localise the impurity state. This phenomenon of self- 
trapping of neutral impurity atoms in quantum degener- 
ate gases has received a significant amount of theoreti- 
cal attention in recent times and our description 
paves the way for an immediate experimental realization 
of the effect in this novel one-dimensional configuration. 
Moreover, our work complements and extends recent the- 
oretical interest in the area of general impurity embed- 
ding in the Tonks-Girardeau regime [I51U7| . 

A low-density gas of N identical bosons trapped in a 
quasi one-dimensional waveguide can be described by the 
Hamiltonian 
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where m is the mass of the particles, z the axial coordi- 
nate, and V{zi) = ^muj 2 z 2 , is the axial trapping poten- 
tial, with W|| the corresponding axial angular frequency 
and ay is the corresponding trap length. The strength of 
the atom-atom contact interactions is given by k, the one- 
dimensional coupling constant, which can be derived by 
a renormalisation procedure from the three-dimensional 
scattering process as K — ih rn ^ D (a± — Ca3£>) _1 [I]- 
Here a± = s/K/mcox is the radial trap width and wj_ the 
radial trapping frequency. The standard s-wave scat- 
tering length is denoted by a^n and C fts 1.4603. For 
the 87 Rb isotope we have a^u ~ 5.3 x 10 _9 m for both 
hyperfine states used in [T|, and m w 1.44 x 10 _25 kg. 
For finite k the Hamiltonian describes an inhomogeneous 
Lieb-Liniger gas and we characterise the strongly interac- 
ing regime by demanding that the contact interaction 
dominates the kinetic energy n 3> h n±u/m, where n\£> 
is the mean linear density of the atoms. In this limit 
the many body problem admits a unique and particu- 
larly elegant solution, as it allows one to replace the 
contact interaction in eq. ([I]) with the nodal bound- 
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can be enforced a priori by Slater determinant factorisa- 
tion, ty F (z 1 ,z 2 ,...,z N ,t) = ^ d ^(^7)=(0,i)^ n ^''^' 
where the ip n are the single-particle eigenstates of the 
non-interacting system. We have adopted the conven- 
tion of labeling the first N single-particle eigenfunc- 
tions with the index n = 0, 1, 2, ... A — 1. This, how- 
ever, leads to a fermionic rather than bosonic exchange 
symmetry, which is corrected by proper symmetrisation, 
A = ]li<i<j<jv sgn(zi-Zj) to give * B = A^ F [5J. From 
this exact solution, the time-dependent single-particle 
density is given by [HI [20] 
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In the experiment by Palzer et al. [T] an array of one- 
dimensional clouds was created with the axis vertically 
aligned. A localized excitation of a small section (near 
the center) of the cloud to an untrapped hyperfine state 
was introduced. While these untrapped atoms fall under 
gravity, they remain strongly confined in the transverse 
direction created by the optical fields. The dynamics of 
the system can therefore be split into the description of 
its two components: a trapped background gas of initially 
iVb atoms, and an impurity wave-packet of iVj atoms, 
with atom densities pt,(z,t) and pi(z,t), respectively. 

Girardeau originally used his famous mapping theorem 
to find the unique solution to the single component many- 
body problem. The solution of the full time-dependent 
two-component many-body problem is challenging and 
only a small number of exact solutions exist for time in- 
dependent mixtures [19]. Here we suggest that a two- 
component system with transport can be treated using a 
phenomenological approach based on the time evolution 
of non-linearly coupled single particle states. The density 
distributions are then built from these states by assum- 
ing single component Tonks- Girardeau gases for both the 
impurity and parent components (see eq. Q). 

At time t = 0, all atoms in the system are assumed to 
be in the trapped TG gas state, the density of which can 
be calculated without approximation using Eq. (g. We 
then simulate the experiment through the application of 
a short, broadband pulse, f(z,t), with a Gaussian spa- 
tial intensity of FWHM a = 2.3/im, for a duration of 
£ P uisc = 200^s. The intensity of this pulse, 7, is chosen 
such that the population transfer corresponds to the ex- 
perimental situation in which, on average, three atoms 
were coupled out [TJ . The integrated densities (norms) of 
the two components vary in time and according to a set 





















.., 








K=0 Jm 






K=1. 35><10 3 


B Jm 




K=1.55*10 3 


Jm 





-20 -10 10 20 
Distance, z[pm] 



500 1000 1500 2000 
Time [ps] 



FIG. 1. Density profiles of the two-component Tonks gas sys- 
tem for different times (a) t = 200 ps, (b) t = 1400 ps and 
(c) t = 2000 ps for N b = 50 87 Rb atoms and cjy = 2tt x 39 
Hz. The grey (lighter) curves show the trapped component, 
while the red (darker) curves show the impurity wave-packet. 
The width of the outcoupling pulse is a — 2.3 pm and the 
inter-component interaction strength is k = 1.35 x 10~ 36 Jm. 
Subfigure (d) shows the center-of-mass trajectory of the im- 
purity atoms for different interaction strengths k. 



of coupled non-linear Schrodinger equations of the form 
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The densities of the impurity and background com- 
ponents are then built according to pi(z,t) — 

J2n=o I'WM)! 2 and pi(z,t) = J2n=o I'M-M)! 2 respec- 
tivley. Here g is the gravitational acceleration and we 
have assumed that the inter-component scattering is in 
mean-field regime. This can be justified by realising that 
the impenetrability of the non-alike atoms in the Tonks 
limit is overcome by the gravitational pull on the im- 
purity, which leads to an increase in the scattering pa- 
rameter fc|ai£)|, where k is the linear atomic momentum 
and am is the atomic scattering length. In this regime 
the collisions between particles in different states can no 
longer be treated as collisions between impenetrable par- 
ticles. The strength of this scattering, n, is calculated 
from the data given in [1] to be n ~ 1.35 x 10~ 36 Jm. In 
the following we will show results from simulating these 
coupled equations for a gas of initially Ab = 50 atoms, 
from which the outcoupling process transfers A; ps 3 into 
the impurity. All parameters are chosen such that a di- 
rect, qualitative comparison to the experiment is possible 

The densities of the components after the outcoupling 
process and during the subsequent dynamics are shown 
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FIG. 2. (a) Atom blockade: shown is the fraction of out- 
coupled atoms remaining within the confines of the trapped 
gas after t — 3000 /is. (b) Density profile taken at t = 2000 /is 
for a gas in the intermediate regime at k — 1.44 x 10 -36 Jm. 

in Fig. [T] Initially, after the coupling pulse is switched 
off (t = 200^s), the impurity component is localised at 
the origin (z = 0[im) and subsequently disperses and dis- 
places as it is dragged through the trapped component 
by gravity (t = 1400/xs and t = 2000/is). The response 
of the background cloud on these timescales is mainly 
determined by the interaction with the falling impurity 
and our results agree qualitatively well with the exper- 
imental findings of Palzer et al. [1]. This provides sup- 
port for our specific phenomenological model, as other 
approaches, for example assuming coherence within the 
impurity, led to a dynamics that agreed very poorly. 

The transport process was studied by analysis of the 
impurity's center-of-mass motion, which is shown in 
Fig. [fjd) for three different values of the interaction 
strength, re. It can be seen that a finite interaction pro- 
duces noticeable deviations from parabolic flight. Firstly, 
the center-of-mass is expelled downwards, on a very short 
time scale. This can be understood by realising that the 
resonance position for the coupling is moved towards neg- 
ative z-values for increasing magnitude of the non-linear 
terms in Eq. ^ (see also discussion below). After this 
the fall of the outcoupled atoms is slowed by the trapped 
cloud as compared to the free-fall case. We find that for 
re = 1.35 x 10 -36 Jm eventually all atoms exit the overlap 
region with the trapped component, however the signif- 
icant slowdown for re = 1.55 x 10~ 36 Jm points towards 
the interesting possibility of blockade within the gas. 

Being able to trap the impurity and leave it embedded 
in the parent cloud is an exciting prospect for a novel 
mixed quantum gas. We therefore calculate the fraction 
of impurity atoms that might be retained within the con- 
fines of the trapped cloud ([— 20, 20]/im) after 3000 [is. 
In each case approximately three atoms are present in 
the impurity wave-packet and Fig. |2ja) shows how the 
fraction remaining within the cloud region varies as a 
function of re. Three regimes can be identified: (i) at in- 
teraction strengths below re w 1.35 x 10 -36 Jm the entire 
impurity wave-packet is able to pass through the trapped 
component and exit the cloud; (ii) in the region between 
re « 1.35 x 10~ 36 Jm and re w 1.5 x 10 -36 Jm the im- 
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FIG. 3. Output yield for a pulse of duration t pu isc = 200/is. 
Panel (a) shows the asymmetry present when the focus of 
the pulse is offset by a distance zq = ±6/mi from the center 
of the cloud. Panel (b) shows how the frequency of the Rabi 
oscillations (at zo = 6/im) are affected for different strength of 
gravity, g c s/g = 0.2, 0.4, 0.6, 0.8, 1 corresponding to the curves 
from top to bottom. The curves are vertically offset for clarity 
by AN = 0.2. Note that the justification for the non-linear 
term in eq. |3} requires a finite gravitational strength. 



purity wave-packet splits into two components of which 
one leaves the cloud and the other remains trapped in 
the center; and (iii) above re m 1.5 x 10 -36 Jm maximum 
self-trapping is achieved, with > 90% of the atomic den- 
sity blockaded inside the trapped gas for long timescales. 
A typical density plot for an interaction strength in the 
intermediate regime (re = 1.44 x 10 - 36 Jm) after 2000/is 
is shown in Fig. [2jb) . Approximately half of the out- 
coupled density is clearly visible to be localised in the 
center of the background gas, while the other half has 
been significantly accelerated by gravity. This part will 
continue to travel out of the cloud. 

Let us in the following consider the effects due to the 
asymmetry induced by the gravitational potential. For 
this we look at the situation in which the centre of the 
outcoupling pulse is located away from the origin of the 
harmonic potential. Fig. [3^a) shows the results for the 
number of outcoupled atoms when the pulse is focused 
around zq — ±6/im. One immediately notes that the out- 
coupling yield at both positions is dramatically reduced 
in comparison with coupling from the center. While this 
can be partly explained by the lower density of the TG 
gas at z = ±6/xm, it alone is not sufficient for the dra- 
matic reduction. In fact, the low efficiency is due to the 
increasingly off-resonant nature of the outcoupling pro- 
cess for larger distances from the centre, which stems 
from the presence of the strong gravitational potential 
in the outcoupled channel. This is also consistent with 
the increased value of the Rabi frequency as can be seen 
from Fig. [3^ a) for zq = 6/im. In fact, carrying out simu- 
lations for a system in which the strength of gravity can 
be decreased (corresponding to the optical lattice being 
rotated from vertical to horizontal), shows in increase in 
Rabi frequencies with decreasing detuning between the 
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channels (see Fig.|3jb)). Note that while during the time 
of the outcoupling process the gravitational pull has only 
a small effect on the position of the outcoupled atoms, 
the dephasing of the output can already be seen in the 
damping of the Rabi oscillations. 

To understand the yield variations observed in 
Fig. [3ja) let us in the following investigate the reaction 
of the system to a change in the width, a, of the outcou- 
pling pulse. In the middle panel of Fig. [4] we show the 
output yield at the end of the coupling pulse (t = 200^s) 
as a function of the position of the focal point for a pulse 
of the experimental FWHM of a = 2.3/im. Focusing on 
or close to the center of the trap ensures the expected 
large efficiencies, and the strong asymmetry observed in 
Fig.[3]ja) becomes visible for increasing values of \z§\. The 
visible fine-structure is due to the existence of Rabi os- 
cillations. The strong fall-off away from the trap centre 
shows that the effect of gravitational detuning is impor- 
tant over the scale of a few microns and therefore can 
be important over the spatial profile of the outcoupling 
pulse. It will, in particular, influence the coherence of the 
coupling process and we show in the left and right panels 
of Fig. g] that an increase of the pulse width, <r, leads to 
a decay of the coherent, high-contrast Rabi oscillations. 
The fact that for outcoupling pulses of the same duration 
the number of oscillations visible for small values of a is 
different for coupling above and below the cloud centre 
can be understood by remembering that the resonance 
point is below the the trap centre. Therefore the detun- 
ing at positive Zq is larger than at negative z , leading to 
higher frequency oscillations. 

As the spot size increases the detuning gradient gives 
rise to different phase components in the outcoupled 
pulse, which leads to the vanishing of the Rabi oscilla- 
tions. This is also the reason for the displacement of the 
oscillations with increasing spot size, as broader spots are 
able to sample areas closer to the resonance (which for 
our parameters lies at zr — — 1.9/im) with lower Rabi 
frequency. As the spot size increases further a pulse 
at zq = — 6/im significantly overlaps with the resonance 
point and its output becomes dominated by a single (co- 
herent) component from the high-density central region 
of the trap. This leads to the increase observed in the 
upper right hand side corner of the left panel of Fig. [4] 
While the spot size for the outcoupling pulse cannot be 
made arbitrarily small, these effects could be experimen- 
tally observed by weakening the longitudinal trapping, 
for example. 

In conclusion, we have presented a model for trans- 
port in strongly correlated quantum gases that is able to 
reproduce key features of recent experiments. Further- 
more, our model predicts several features which could 
be explored further experimentally. We have brought 
strong evidences that even in the presence of gravity a 
strong interaction between an out-coupled and a residing 
component in a Tonks gas can lead to a self localization 
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FIG. 4. The effect of pulse width a on the atom output yield 
from a TG gas below and above the centre of the trap. Left: 
atom yield for a pulse focused at zq — — 6/im. Center: Output 
yield (at the end of the pulse) as a function of focal point for 
a — 2.3/im. Right: yield for focus at zo = 6/im. Well-defined 
Rabi oscillations are clearly visible for small a for for both 
situations, above (right) and below (left) the trap centre. 

of the out-coupled component. Although atom blockade 
was not observed in the experiment of Palzer et al., the 
regime could be reached by tuning a Feshbach resonance 
between the hyperfine states. In addition our work has 
shed light on interesting aspects of the out coupling pro- 
cess which are unique to the configuration under consid- 
eration. Using our model to investigate other condensed 
matter phenomena in these correlated systems, such as 
spin-charge separation and analogues of Cherenkov radi- 
ation, are subjects of ongoing research. 
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